Stress-free Spatial Anisotropy in Phase-Ordering 
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We find spatial anisotropy in the asymptotic correlations of two-dimensional Ising models under 
non-equilibrium phase-ordering. Anisotropy is seen for critical and off-critical quenches and both 
conserved and non-conserved dynamics. We argue that spatial anisotropy is generic for scalar 
systems (including Potts models) with an anisotropic surface tension. Correlation functions will 
not be universal in these systems since anisotropy will depend on, e.g., temperature, microscopic 
interactions and dynamics, disorder, and frustration. 
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Most theoretical, numerical, and experimental treat- 
ments of non-equilibrium phase-ordering kinetics assume 
that asymptotic correlations are spatially isotropic in 
stress-free systems quenched from disordered initial con- 
ditions into the ordered phase . Numerical measure- 
ments in simple Ising models have supported this as- 
sumption! B|4j, despite several qualitative reports to the 
contrary [gl]. It has not been clear whether Potts models 
and frustrated 3d Ising models where anisotropy 
effects have been measured, are typical or are special 
cases. Indeed, without a demonstration that anisotropy 
is expected asymptotically it has been possible to discount 
measured anisotropy as being a transient effect (see, e.g., 

§)• . _ 

Naively, it is reasonable to expect anisotropic correla- 
tions in phase-ordering. The 2d Ising model below Tc, 
for example, is spatially anisotropic at arbitrarily large 
scales in equilibrium correlations |^ , in the resulting sur- 
face tension , and hence in the interface dynamics . 
In fact, we demonstrate below that spatial anisotropy is 
generic in scalar systems for quenches into an ordered 
phase with an anisotropic surface tension. Isotropic the- 
ories are inadequate for these systems. 

For a coarse-grained scalar order parameter, (p, in the 
continuum limit, we define an effective free energy in mo- 
mentum space = / (fik[{k'^ + £'k)(/'k0-k + Vfe], 
where [{</>}] is the fourier transform of local poten- 
tial terms and the differential operator includes 
anisotropic higher-order gradients. We define an angle 
dependent free-energy density e(n) by restricting the in- 
tegral in F to momenta in a given direction n and aver- 
aging (denoted by angle brackets) over the random initial 
conditions. To do the restricted momentum integral in 
{F) , at late enough times for domain walls to be well de- 
fined, we use the anisotropic Porod law for the structure 
factor S'(k) = {4>k(t>^k) in general dimension d (general- 
izing PI): 



(1) 



orientation n = k/fc. A{t) ~ L^^ is the average area 
density of domain wall and P(n) is the angle distribution 
function of domain wall orientation [with P(n) — P(— n) 
and / dnP(n) 1]. We find 



e(n) = Aa(n)P(n), 



(2) 



where the leading k^ and potential terms make isotropic 
contributions to (T(n), while makes anisotropic con- 
tributions both directly and through ^(n)~^. From Eqn. 
(^), we identify a{n) as the effective angle dependent 
surface-tension for domain walls with normal direction 
n. 

Now consider a dissipative quench to T = 0, with 
no thermal noise. The dynamics will be given by 
0k — — r(k)(5P/(50_k, where the dot indicates a time- 
derivative, r = const, corresponds to non-conserved dy- 
namics and r = T2k^ to conserved dynamics p^ . The 
time rate of change of the energy-density is then simply 



e(n) = J dkk'^-\SF/S(Pk^k) 
= - I dkk''-^T-\^^^_^) 



which holds for ^ |k| ^ ^(n)~^, where L{t) is the 
characteristic growing length-scale of the system and ^ (n) 
is the domain wall width for domain walls with normal 



(3) 



where k = kn. We have used the dynamics of to replace 
the functional derivative of the free-energy so that the 
resulting expression for e(n) has no explicit dependence 
on F[{^}]. 

Under the assumption of isotropic correlations, we 
generically obtain different anisotropics in Eqns. (|^) 
and (^) at arbitrarily late times — a contradiction since 
e = dtc. Apart from possible anisotropics in the correla- 
tions, the anisotropy of e(n) in Eqn. (H) is determined 
by the statics, through o'(n). On the other hand the 
anisotropy of e(n) in Eqn. (^) is determined by the dy- 
namics, through r(k), in addition to possible contribu- 
tions by the effective UV cutoff ^(n)^^. Since a and F 
are independent, the anisotropics will not generally be 
equal unless anisotropic correlations make up the differ- 
ence. For the general case, anisotropy in cr(n) implies 
anisotropy in S'(k) even in the scaling limit J15[. 
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The renormaUzation-group (RG) approach to phase- 
ordering [|l6| is easily generalized to include anisotropy. 
The only change is to note that any anisotropy of either 
-F[{0}] or r(k) will be renormalized by microscopic de- 
tails. (An illustration of this renormalization is the tem- 
perature dependence of the effective surface tension Q.) 
The demonstration that thermal noise will be asymptot- 
ically irrelevant for quenches to below Tc will still ap- 
ply, with the caveat that the effective T = dynamics 
will include the surface-tension at the quench tempera- 
ture. We then apply our above argument that predicts 
anisotropy with noise free dynamics. As a result, we 
expect anisotropy for all scalar quenches below Tc WA - 
Anisotropy may be renormalized by temperature, disor- 
der, the details of the local interactions in the system 
(such as frustration), and even by global conservation 
laws that are "irrelevant" in terms of growth laws. 
Anisotropy may also depend on the details of the micro- 
scopic dynamics. 

In principle, we could try to choose the anisotropy of 
r(k) to allow isotropic correlations despite an anisotropic 
cr(n). This fine-tuning of the dynamics (F) with re- 
spect to the statics [a) will not generically occur. The 
RG approach shows that r(k) will only be renor- 
malized analytically, i.e. anisotropy will only enter at 
0{k'^) and higher. For conserved dynamics, these contri- 
butions are subdominant in Eqn. ^ since the integral 
converges in the UV |l^] . Thus neither the anisotropy of 
F(k) nor the anisotropy of the core scale ^(n) will affect 
e(n) through Eqn. (^). For conserved dynamics, even 
fine tuning of F cannot eliminate anisotropic correlations. 
With non-conserved dynamics, the UV regime dominates 
the energy-dissipation integral (|^) and both F(k) and 
^(n) make anisotropic contributions to e(n). In principle 
the anisotropy of F(k) could be renormalized to compen- 
sate cr(n) — allowing isotropic correlations within our 
argument. However, we find anisotropy numerically for 
T > even for non-conserved dynamics. We deduce that 
fine tuning of the kinetic prefactor does not occur. 

For a system defined on a lattice, anisotropies will be 
present in the surface tension at T = 0, because lattice 
interactions are not rotationally invariant. Our argument 
then implies anisotropic correlations in the scaling limit. 
Anisotropy is only implied in some of the correlation 
functions and not necessarily in two-point correlations. 
In practice we find, in all of our numerical studies, sig- 
nificant anisotropy effects in the two-point correlations. 

For the remainder of this letter, we explore 2d Ising 
models with nearest-neighbor interactions on a square 
lattice under quenches from random initial conditions. 
We find anisotropic correlations quite generally — for a 
variety of temperatures below Tq, of initial magnetiza- 
tions, and for all of globally conserved, non-conserved, 
and locally conserved dynamics. These anisotropies do 
not decrease at late times, as would be expected for tran- 
sient effects introduced by the dynamics at earlier times. 




FIG. 1. A 512^ region of a globally conserved T /Tc = 0.2 
quench with {</>) = 0.4 (t — 513 MCS). Lattice directions, 
here and in the next figure, are vertical and horizontal. 

We measure the normalized correlations C(r, t) = 
m^)m) - (0)^)/((</'(O+)0(O)) - (</))2), which ranges 
from 1 at short distances to at infinity The 
anisotropic length-scale L(n, t) of a system is defined 
by the scale in direction n at which C = 0.5 for non- 
conserved dynamics, and by the first zero of C for con- 
served dynamics. We scale correlations in all directions 
by the length-scale in the diagonal direction. A natural 
measure of anisotropy is x = {Lmax/Lmtn - 1)/(V2 - 1), 
where Lmax is the maximum length-scale at a given time, 
Lmin is the minimum, and the normalization makes x run 
from (circle) to 1 (square) for convex contours of C(r). 

We first consider off-critical quenches with a global 
conservation law to prevent the magnetization from sat- 
urating. We couple the system to a Creutz spin reservoir 
of size 2 each randomly chosen spin is updated by 
a Metropolis algorithm, subject to an additional micro- 
canonical constraint that any spin change (±2) fits in 
the spin reservoir. We study size 1024^ systems with 
((/>) = 0.4. A snapshot from a quench to T = 0.2Tc in 
Fig. |l| illustrates the strong anisotropy even above the 
roughening transition p9[ | . We show some contour plots 
of the scaled correlations in Fig. ^. It is clear that the 
anisotropies are not limited to the small r regime. The 
anisotropy is increasing at late times (see inset of Fig. |^). 
In the same regime, the spherically averaged correlations 
scale well. The latest anisotropies, &t t — 2049 MCS, 
before finite-size effects entered were x = 0-45 (T = 0), 
0.38 (T = 0.2Tc), and 0.12 (T = 0.4rc). [Statistical error 
bars, with at least 30 samples in each case, are less than 
±0.001.] 
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FIG. 2. Anisotropic contours of scaled correlations 
C(r/L) = 0.9, 0.8, 0.7, 0.6, and 0.5 (from the center) for 
an off-critical quench to T = with {(p) — 0.4 and globally 
conserved Creutz dynamics. The times are t = 513 MCS 
(dotted), 1025 (dashed), and 2049 (solid). Also shown, scaled 
by 1.4 for clarity, are the C — 0.5 contours of quenches to 
T/Tc = 0, 0.2, 0.4, and 0.9 (solid, dashed, dot-dashed, and 
dotted lines, respectively) at t = 1025 MCS. The length-scale 
L is such that, along the diagonal direction, C{L) — 0.5. 

We also studied non-conserved critical quenches. With 
heat-bath dynamics and a sublattice update, late times 
in large lattices could be reached. Even so, the asym- 
metry remained small. For lattices of size 2048^, and 
a quench to T = 0, we show dC/dx vs. x {x = r/L), 
along the diagonal and lattice directions, in Fig. || [pO[ . 
We find X = 0.03 at the latest time {t = 4097 MCS, 74 
samples, statistical error ±0.0003). The constant but ori- 
entation dependent domain wall width, ^(n), evident in 
the sharp downturn of dC / dx near a; = in the diagonal 
correlations, increases x by 0{1/L). This is significant 
for small anisotropics and early times. Directly subtract- 
ing this 0(1/L) contribution leads to x slowly increasing 
with time (inset of Fig. ^), with a corrected latest value 
X = 0.02. 

We have also simulated conserved 2d Ising systems 
with nearest-neighbor Kawasaki exchange dynamics and 
a Metropolis update. At low enough temperatures for 
the anisotropy of cr(n) to be visible, the activated dynam- 
ics slows the simulations considerably. We explored size 
256^ systems, with (0) = 0.4 and T = 0.4rc, up to times 
t — 10^ MCS (10 samples). The length scales achieved 
(L 12) are so small that x ~ within numerical accu- 
racy, so in Fig. ^ we plot C(r) against the energy-energy 
correlation function C£;(r) = (£;(r)i5(0))/(ii^)^-l, where 
E{y) is the number of broken bonds at site r minus the 
equilibrium bulk average. This shows a significant and 
increasing difference between correlations in the lattice 
and diagonal directions. 
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FIG. 3. dC /dx vs scaled distances, x, for a non-conserved 
critical quench to T = 0. The plusses indicate correlations in 
the lattice axis direction, while the triangles indicate corre- 
lations along lattice diagonals (at t = 4097 MCS). Solid and 
dashed lines indicate corresponding correlations at t = 2049 
and 1025 MCS respectively. In the inset is the anisotropy 
measure x vs. t for T /Tc = 0, 0.2, and 0.4 globally-conserved 
quenches from Fig. 2 (crosses, dotted lines, and diamonds, 
respectively). At the bottom of the inset are the bare and 
corrected x (solid and dot-dashed lines, respectively) for the 
non-conserved quench of the main figure. 
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FIG. 4. Two-point correlations vs energy-density en- 
ergy-density correlations for a conserved quench to 
T/Tc = 0.4 with (0) = 0.4 (size 256^). The correlations are 
spherically averaged (stars and circles), along the axis (plusses 
and squares), and along the diagonal (crosses and triangles). 
For clarity, solid lines have been used after the first zero of C. 
Times are 2.6 x lO'^ and 1.0 x 10*^ MCS, respectively. Dotted 
lines show data from size 128^ systems (36 samples) at the 
earlier time. 

In summary of our numerical results, we find 
anisotropic correlations in various quenched 2d Ising 
models. Anisotropy increases with decreasing temper- 
ature and for increasing net magnetization. Anisotropy 
effects are always slowly increasing at the latest times 
of our simulations. In all of our simulations, the spher- 
ically averaged correlations scale reasonably well while 
the anisotropy is still evolving. To study the non-zero 
asymptotic anisotropy in these systems, some sort of ac- 
celeration method is needed (see, e.g., [|j) — though in 
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general the anisotropy will depend on the numerical al- 
gorithm used. 

In disordered HJI^ and frustrated models, it has 
been argued that scaling functions will be "universal" 
— identical to those of the Ising systems that we have 
studied in this letter. [While logarithmic growth is seen, 
it is thought to come from L dependence in the kinetic 
prefactor F — so that scaled correlations are unaffected.] 
However, anisotropy should be renormalized by disorder 
and frustration, so we do not expect scaling function uni- 
versality to hold in these systems. In frustrated 2d and 3d 
Ising models, the fairly large anisotropy seen numerically 
should remain at arbitrarily late-times jl^. Hope- 
fully, in experimental random- field systems (e.g. [^ ), 
anisotropy can be measured directly. 

We can generalize our argument around Eqns. (|l|)-(||) 
to systems with other types of singular defects. Using 
vector 0{n) order-parameters, and a generalized Porod's 
law S'(k) = D(k/A:)L-"A:-('*+") §, we find that the 
anisotropic contribution to the energy density is asymp- 
totically negligible for systems without domain walls . 
However, other systems with dissipative dynamics in 
which domain walls dominate the asymptotic energetics 
will be anisotropic if the surface tension is anisotropic, 
e.g. Potts models (see 0). 

The growth laws of the characteristic length scale L{t) 
will remain independent of any anisotropics present, as 
long as dynamical scaling is maintained. This follows 
from the Energy-Scaling approach [ p^ since anisotropy 
docs not change the scaling properties of the energy 
or the rate of energy dissipation. We would be sur- 
prised if the anisotropy affected the dynamical scaling 
of the correlations (see however [^), though the scaling 
regime seems to be pushed to much later times as the 
anisotropy slowly develops. It remains an open question 
whether non-zero anisotropics have implications beyond 
the scaled correlations, such as in autocorrelation expo- 
nents. 

In practice, isotropic theories have worked fairly well 
for spherically averaged correlations. Certainly, lattices, 
interactions, and dynamics can be chosen to minimize 
anisotropics. This would be desirable, for instance, in 
lattice simulations of isotropic fluid or polymer systems. 
However, the language of an isotropic zero-temperature 
phase-ordering fixed point is inappropriate for a scalar 
system with an anisotropic surface tension. 

In summary, we expect anisotropy for any scalar lat- 
tice system quenched to below Tc, including disordered 
and/or frustrated systems. We expect the anisotropy to 
depend on the details of the system. Scaled correlation 
functions of anisotropic systems will not be universal. 

I thank J. Cardy for discussions. This work was sup- 
ported by EPSRC grant GR/J78044. 
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